Topologically trivial zero bias conductance peak in semiconductor Majorana wires 

from boundary effects 
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We show that a topologically trivial zero bias conductance peak is produced in semiconductor- 
superconductor hybrid structures due to a suppressed superconducting pair potential and/or an 
excess Zeeman field at the ends of the heterostructure, both of which can occur in experiments. The 
zero bias peak (ZBP) (a) appears above a threshold parallel bulk Zeeman field, (b) is stable for a 
range of bulk field before splitting, (c) disappears with rotation of the bulk Zeeman field, and, (d) 
is robust to weak disorder fluctuations. The topologically trivial ZBPs are also expected to produce 
splitting oscillations with the applied field similar to those from Majorana fermions. Because of 
such strong similarity with the phenomenology expected from Majorana fermions we find that the 
only unambiguous way to distinguish these trivial ZBPs (of height 4e 2 //i) from those arising from 
Majorana fermions (of height 2e 2 /h) is by comparing the (zero temperature) peak height and/or 
through an interference experiment. 



The search for zero-energy Majorana fermions [IH1] 
in solid state systems has received tremendous interest 
[SHE] in recent theoretical and experimental studies. A 
practical route to engineer such a system in the labora- 
tory has been suggested in Refs. [5HT2"]. following earlier 
similar proposals in topological insulators [13] and cold 
fermion systems [14j[15]. The proposed system consists of 
a spin-orbit coupled semiconductor thin film or nanowire 
proximity coupled to a s-wave superconductor and in 
the presence of a suitably directed Zeeman spin split- 
ting. The applied Zeeman field drives the engineered hy- 
brid nanowire (semiconductor Majorana wire) through a 
topological quantum phase transition (TQPT) to a topo- 
logically nontrivial superconducting phase with localized 
zero energy Majorana bound states (MBSs) when the 
Zeeman splitting Y exceeds a critical value T c . The Ma- 
jorana wire is in a topologically trivial superconducting 
phase with no defect-localized MBS when the Zeeman 
splitting satisfies Y < Y c . The zero energy MBS local- 
ized at the ends of the nanowire (for Y > Y c ) has been 
proposed to be observable by a zero bias conductance 
peak in charge current through the ends of the semicon- 
ductor wire [lOl [P6H27] . Here we show that such a charge 
current zero bias peak (ZBP) can exist even in the topo- 
logically trivial phase (r < Y c ) due to excess Zeeman field 
and/or suppressed superconducting pair potential at the 
wire ends, both of which can occur in experiments. Our 
calculations establish that a (non-quantized) ZBP from 
tunneling into the ends of a spin-orbit coupled wire is 
unable to produce an unambiguous signature of MBS, 
even after taking into account all the proposed experi- 
mental checks for Majorana fermions that have appeared 
in recent theoretical and experimental works [551 f2"?H2"5] . 

The critical Zeeman field Y c (= yj A 2 + /i 2 ) in the 
quasi- ID heterostructure depends on the value of the su- 
perconducting pair potential A proximity induced in the 
semiconductor and the chemical potential /i of the semi- 



conductor measured from the top-most occupied confine- 
ment energy band (by "band" , we mean a pair of spin- 
split sub-bands, the band themselves being separated by 
the energy gaps due to lateral confinement). The value 
of r c (for A — > 0) corresponds to the value above which 
(r > /i) an odd number of confinement induced sub- 
bands are occupied in the semiconductor. For an even 
number of sub-band occupancy the system is in the topo- 
logically trivial phase with no MBS localized at the ends 
of the nanowire. 

The typical experimental nanowires, however, are more 
likely to have an even number of sub-band occupancy 
than odd. This is because the confinement energy gap 
Ec is typically much larger than the Zeeman energy gap 
P, and, consequently, the chemical potential is expected 
to lie in the confinement gap between two pairs of spin- 
split sub-bands rather than in the Zeeman gap between 
two sub-bands in the same band. By a simple estimate 
(within the non- interacting electron model), assuming 
that the ID density of states cx 1/VE, the ratio of the 
number of samples with an even sub-band occupancy to 
that with an odd sub-band occupancy is ~ \jEcjY ~ 6. 
Here we have used typical values for Ec ^6 — 7 meV 
and the applied Y ~ A ~ 150- 20(hieV from Ref. [25]. It 
follows that, unless /i can be tuned by an externally ap- 
plied gate potential (which is hard because of the strong 
coupling to a superconductor) the chemical potential in 
the experimental nanowires may cross both sub-bands of 
the top occupied band and the system may be topolog- 
ically trivial. Here we show that, even in this case, a 
robust but topologically trivial ZBP can still appear in 
tunneling experiments due to excess local Zeeman field 
and/or suppressed superconducting pair potential at the 
wire ends. The excess Zeeman field at the wire ends can 
be due to the Meissner effect (the edge magnetic field 
larger than that in the bulk), or due to the presence of 
a magnetic impurity at the nanowire-lead interface, and 
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the suppressed A at the ends of the wire can be a con- 
sequence of the (spatially varying) proximity effect near 
the boundary of the bulk superconductor. Thus, both 
effects, capable of producing a ZBP even in the topologi- 
cal^ trivial phase, can occur in the experimental systems. 

About the origin of the topologically trivial peaks in 
our work we speculate on the following scenario: As dis- 
cussed in Ref. [30] (see also Ref. [55]) the zero energy 
bound states from the two different sub-bands (in the 
top occupied band) at the same end of the nanowire could 
contribute to a total ZBP height of 4e 2 //i provided their 
coupling and the resultant energy splitting could be sup- 
pressed. Ref. [3U] proposed that a large uniform bulk 
Zeeman field, coupled with a smooth barrier potential at 
the lead-nanowire interface, can suppress the inter-sub- 
band s-wave pairing [31] , resulting in a near-zero-energy 
peak even in the topologically trivial phase. We specu- 
late that the coupling between the two zero energy states 
can be suppressed even when a perturbation is applied 
only near the heterostructure ends because the zero en- 
ergy states from the (uncoupled) sub-bands are localized 
there within a length scale ~ £, the coherence length. 
The local excess Zeeman field and/or suppressed pair 
potential at the nanowire ends can however be a nat- 
ural consequence of the experimental geometry and thus 
present an alternative explanation of the (unquantizcd) 
ZBPs seen in the recent experiments on semiconductor- 
superconductor heterostructures [2"5H2"7] 

The trivial ZBP appears in our calculations entirely in 
the topologically trivial phase, that is, with only an even 
number of confinement induced sub-bands occupied. Ac- 
cording to the scenario presented above, even though the 
ZBP (of height Ae 2 /h) can be viewed as the result of 
resonant local Andreev reflection from two (or an even 
number of) weakly coupled MBSs, it is important to men- 
tion that no TQPT is necessary for the near-zero-energy 
states to appear. In other words, an even number of Ma- 
jorana fermions, even if they are only weakly coupled, 
can be thought of as a conventional Dirac fermion, and 
do not have non-Abelian statistics. It is in this sense 
that we call the ZBPs observed with local perturbations 
in our work topologically trivial ZBPs. 

We consider a one-dimensional multichannel Rashba 
spin-orbit coupled semiconducting nanowire (for exam- 
ple, InSb) oriented along the x-direction, and in prox- 
imity to a s-wave superconductor (for example, NbTiN). 
The full hybrid structure is modeled by a discrete N x W 
square-lattice tight-binding model [32l [33] 

Hnw = Hke + H so + Hz + H sc , (1) 

H K E = ^(<4+d,a C r,a + H - C + tt4,a c r, a ~ 

r,d,a 

Hso = ^ -iaRcl +dt0t z.(a a ^ X d)c r)/ 9 + H.C., 

r,d,a,/3 

Hz = c i,a( B - a )^^ H sc = - Ac tA-.i + H ' C " 

r,a,/3 r 

where r G {l,m} denotes lattice sites with the index I 



along the length, I — 1,2..N, and the index m along 
the width of the nanowire, m = 1,2..W, (W counts the 
number of tight-binding chains) d € {e x , e y } is a unit vec- 
tor connecting nearest-neighboring sites, (a, f3) € {\, 1} 
are spin indices, and er is spin-1/2 Pauli matrix vec- 
tor. Here c\ a is a creation operator of an electron with 
spin a at lattice site r of the nanowire. The hopping 
strength t of the electrons is related to the band mass 
m* = h 2 /(2ta 2 ) where a denotes lattice spacing in the 
tight-binding model. The applied magnetic field B g 
{B x , By, B z } opens a Zeeman splitting T — (gfig/2)B in 
the sub-bands, and A is the proximity induced supercon- 
ducting pair potential in the semiconductor. In the fol- 
lowing we assume gfis/2 — 1 and identify T with B in the 
reduced units. In recent experiments [2"oj|"2T] B along the 
wire axis (a;-axis) is increased to observe a ZBP, and then 
B is tilted from the wire axis to remove the ZBP, both 
observations consistent with the Majorana origin of the 
ZBP at the wire ends [TU]. In Eq. fl = 2{^ — t) and an 
energy shift — jl/2 has been added for the local Majorana 
transformation used in our transport calculations. We 
couple the semiconductor-superconductor hybrid struc- 
ture in Eq. ([lj to two metallic leads at the two ends. The 
metallic leads are modeled by free electron tight bind- 
ing chains. Each first {l,m} and last sites {N, m} with 
rn = 1,2..W of all the transverse tight-binding chains of 
the semiconductor nanowire are coupled to a semi-infinite 
free-electron tight binding chain which forms the metal- 
lic leads. All the tight binding chains in the left lead are 
kept at chemical potential fi^ and temperature T^, and 
those in the right lead are at chemical potential and 
temperature Tr. The Hamiltonian H^ (p — L,R) below 
describe the metallic lead Hamiltonians and H^ the cor- 
responding tunnel couplings between the leads and the 
nanowire, 

W oo 

a,m— 1 k=\ 
W 

— ~7p J]] ( c m,l,a c m,l p ,a + c m ,l p ,a ^rn,\,u)- (2) 
a,m— 1 

where l\, = 1 and = N. Here c^J k a is an elec- 
tron creation operator on the pth lead. The strength 
of tunnel coupling between the normal leads and the 
nanowire is 7^. It controls the width and the height of the 
ZBP. Here we apply the quantum Langevin equations and 
Green's function (LEGF) approach [34 38 to calculate 
the current-voltage (I-V) characteristics and the corre- 
sponding differential conductance dl/dV across the mul- 
tichannel spin-orbit coupled semiconductor nanowires. 
To mimic the recent experimental conditions [25j . we set 
Hr ; = and drive /il from — eV to eV . We fix the tem- 
perature of the two leads to be the same, Tl =Tr = 0. 

An applied magnetic field along x splits the two sub- 
bands of a single tight binding chain, and the MBSs ap- 
pear when only one sub-band is occupied, i.e., the chem- 
ical potential lies in the gap between the two split sub- 
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FIG. 1. Zero-temperature differential conductance dl/dV 
vs. the applied voltage V for different numbers of transverse 
chains W for /i = 0. Everywhere B x — 0.8, B y = B z = 0.0 
and |A| = 0.6, so B x > B x<c = |A|. Number of chains 
W in the panels is as follows: (a) W = 1, (b) W — 2, (c) 
IF = 3, and (d) W = 4. Since /tt = in Eq. (fit (with periodic 
boundary conditions) corresponds to an equal number of sub- 
bands below and above E = fi, only W = 1, 3 has a ZBP for 
B x > B XiC = |A|. For W = 2,4, B x > |A| produces no ZBP 
because only an even number of sub-bands are occupied. 



bands. The critical magnetic field along x direction is 
given by B c = ^A 2 + (j, 2 , the MBSs appear for B x > B c . 
The above formula for B c holds even for multichannel 
nanowires, however the chemical potential /i is deter- 
mined from the bottom of the topmost occupied band. 
In our Hamiltonian (Eq. ([!])) /i = lies in the middle 
of all the sub-bands, and thus an odd number of sub- 
bands are occupied (for B x > /i) when the number of 
chains W is odd, and an even number of sub-bands are 
occupied (for B x > /i) when W is even. Therefore, a 
TQPT from a trivial superconducting phase to a topo- 
logical superconductor occurs by increasing B x above B c 
for an odd W but not for an even W (for fj, = 0) . We use 
Fig.fijto illustrate this: there is a ZBP from the MBS for 
W = 1,3 and no ZBP in the topologically trivial phase 
for W = 2,4. We use N = 40, j p = t = 1, i p = 0.2 
(j) = L, R), n = 0, a R = 0.2, B y = B z = and |A| = 0.6 
in all the figures in this paper, if not explicitly stated 
otherwise. In physical units these parameters for a 2 \xm 
long wire correspond to a Rashba spin-orbit coupling (in 
the continuum model) a = a^a — 0.1 eVA, pair poten- 
tial A = 0.3 meV and Zeeman splitting T = 0.4 meV. 
We have confirmed that the results in this paper are ro- 
bust to variations in the values of these parameters in 
particular to the Rashba spin-orbit coupling parameter 
a (Fig.^). In all the figures we quote the parameters in 
units of t = 1. 

Now we turn on a local perturbation to the magnetic 
field at the first sites {l,m} (m = 1,2.. IF) of all the 
transverse chains of the semiconductor nanowire. We 
find that a ZBP appears in the differential conductance of 





FIG. 2. Emergence of ZBP in the topologically trivial super- 
conducting phase (even band occupancy) due to excess Zee- 
man field and/or reduced pair potential at the wire ends. In 
all the panels, W = 2, A = 0.6, and [i = 0. In panels (a,b,d) 
an — 0.2, B x = 0.8. Note that for these values of B x and 
A there is no ZBP for W = 2 (even band occupancy) in the 
absence of a local perturbation (Fig. fija). The local end per- 
turbations producing the topologically trivial ZBPs are as fol- 
lows: (a) 8B x {l,m} = 0.45, (b) S^,f{l,m} — 1, 5/i^{l,m} = 
0.2, and (d) <5A{l,m} = —0.35, where m = 1,2. In panel 
(c) we show, for illustrative purposes, the emergence of the 
trivial ZBP even for a larger value of the spin-orbit coupling, 
an = 0.5 and B x = 1. For local end perturbations we have 
used, 8B x {l,m} = 0.7, 8B X {2, m} — 0.5, that is, the end 
perturbation is applied on two sites and gradually peters of 
towards the bulk. 



the nanowire even in the topologically trivial phase (even 
number of occupied sub-bands for A — > 0) when the mag- 
nitude of B x locally at {l,m} is larger than that in the 
bulk (Fig. 2a). The ZBP also appears in the topologi- 
cally trivial phase when there is a local B z field at the 
first sites (in addition to bulk B x ) of the same order of 
B x (Fig. [2]d). Note that a local B z field is equivalent 
to a local shift of the chemical potential of the electrons 
with one spin with respect to the other. A topologically 
trivial ZBP appears also for a local suppression in the 
superconducting pair potential A at the first sites of the 
transverse chains (Fig. [2H). The heights of the trivial 
ZBPs are of the order 4e / h (Fig. [2| , which is double the 
height produced by the MBSs (Fig.[l]). 

Applying the local perturbation at multiple boundary 
sites on all the chains (Fig. [2J;) we find that the topo- 
logically trivial ZBP appears in a larger range of the lo- 
cal perturbations than when the perturbation is applied 
at a single boundary site (as in Fig. J2^i{2]d{2|1). This is 
also more natural experimentally since the excess Zee- 
man field and/or the suppressed pair potential at the 
wire ends are expected to gradually shift to the bulk val- 
ues in the interior of the nanowire. In Fig.[3^j3]3 we show 
that the trivial ZBPs due to a fixed set of local perturba- 
tions in the Zeeman field (applied at the two end sites) 
appear only above a threshold uniform Zeeman field in 
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FIG. 3. Dependence of the topologically trivial ZBP on the 
uniform Zeeman field B x for a fixed perturbation at two end 
sites on all chains (SB x {l,m} = 0.32 and 8B x {2,m} = 0.2). 
All panels have W = 2 (two transverse chains) and A = 0.6. 
The ZBP is split in panel (a) (B x = 0.7), is fully developed 
in panel (b) (B x = 0.8), and stays in panel (c) (B x — 1.0), 
before splitting again with larger B x (not shown). Panel (d) 
(B x — 0.4, By = 0.4) shows the disappearance of the ZBP 
with rotation of the uniform field in the plane formed by the 
wire and the effective spin-orbit coupling {{x — y) plane). 
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FIG. 4. Robustness of the topologically trivial ZBP to mod- 
erate disorder fluctuations and the role of barrier potential 
(tunnel coupling with the metallic lead). In all panels W = 2, 
B x = 0.8, and /i-|-{l,m} = 1, ui{l,m} = 0.2. (a) Random- 
ness in /i is chosen from a uniform distribution between and 
0.1, 7p = 0.2 and (b) clean wire (u = 0), but a stronger cou- 
pling with the lead, y'p = 0.3 (p = L,R). These figures should 
be compared with that in Fig. [2Jd. 



the bulk. In Fig. [3]: we show that the ZBP is robust to 
a further increase in the bulk field, before splitting into 
two peaks above a higher threshold (not shown). This 
sort of dependence of the topologically trivial ZBP on 
the uniform bulk field is consistent with the recent ex- 
periments [2"5ll27j . In Fig. |3ji we show that the trivial 
ZBP disappears as the bulk Zeeman field is rotated by 
7r/4 from the axis of the wire in the plane formed by the 
wire axis and the direction of the effective spin-orbit field 
(the (x — y) plane) . This behavior is also fully consistent 
with the recent experiments |25[ I27j , which is taken as 
an argument in favor of the Majorana origin of the ZBP 
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We have checked the robustness of the topologically 
trivial ZBP to a moderate and random impurity distri- 
bution chosen independently for each site, identically dis- 
tributed according to a uniform distribution. We find 
that the trivial ZBP survives even in the presence of 
moderate disorder fluctuations as shown in Fig. [4^,. With 
respect to the clean wire, the height (width) of the triv- 
ial ZBP decreases (increases) with disorder compared to 
that without disorder in Fig.[2}3. However the trivial ZBP 
splits into two peaks for further increase of the strength 
or randomness of the impurity potential. For the depen- 
dence of the ZBP on the barrier potential at the lead- wire 
interface [SHUSH], which, in our calculations, controls the 
coupling 7^ (see Eq. (2)) of the nanowire to the metal- 
lic leads, we find that the height and width of the ZBP 
are somewhat susceptible to the barrier potential. For 
weaker coupling to the lead, the height and width of the 
peak are reduced, while they increase with increasing the 
coupling to the lead as shown in Fig. [4Jd. For a weaker 
coupling to the leads (higher barrier potential) we find 
that the range of local perturbations in which the trivial 
ZBP is formed is reduced compared to the case with a 
stronger lead coupling (lower or smoother barrier poten- 
tial). For the so-called gap-closing signature [2"51 |2"5I I3T)] 
before the appearance of the ZBP, note that no gap needs 
to close before the appearance of the topologically trivial 
ZBP because the system remains in the same topologi- 
cally trivial superconducting phase. This is one impor- 
tant difference between the Majorana origin of the ZBP 
and the topologically trivial ZBP as found here. Even if 
there is no bulk gap closing involved (no TQPT), for a 
linear dependence of the end perturbations on the bulk 
parameters, the spectral weight should continuously shift 
from higher subgap energies to zero energy as the uniform 
bulk field crosses the threshold (Fig. [3]j, [3)3)). This may 
appear as a gap closing signature similar to the calcula- 
tions in Ref. [28 40 . However, for a non-linear depen- 
dence of the local perturbations on the bulk parameters, 
the downward shift of the spectral weight may not be 
gradual with the uniform field, in which case any gap 
closing signature may be suppressed. 

In summary, we investigate the effects of an excess Zee- 
man field and/or suppressed pair potential at the ends 
of the semiconductor Majorana wire on the tunneling 
conductance through the wire ends. Both effects may 
be present in the experiments and we show that both 
can produce stable ZBPs above a threshold uniform bulk 
Zeeman field, even when the system is topologically triv- 
ial (sub-band occupancy even). The trivial ZBPs, since 
they arise from conventional BdG states, are expected to 
show splitting oscillations with the bulk field [55] similar 
to the Majorana fermions. Given such strong similari- 
ties with the ZBPs observable from the MBSs, the only 
unambiguous way to distinguish the trivial ZBPs found 
here (of height 4e 2 /h) from those arising from Majorana 
fermions (of height 2e 2 /h) is by comparing the (zero tem- 
perature) peak height and/or by an interference experi- 
ment 1411. 
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